Nonlinearity and disorder: Classification and stability of nonlinear impurity modes 
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We study the effects produced by competition of two physical mechanisms of energy localization 
in inhomogeneous nonlinear systems. As an example, we analyze spatially localized modes supported 
by a nonlinear impurity in the generalized nonlinear Schrodinger equation and describe three types 
of nonlinear impurity modes — one- and two-hump symmetric localized modes and asymmetric 
localized modes — for both focusing and defocusing nonlinearity and two different (attractive or 
repulsive) types of impurity. We obtain an analytical stability criterion for the nonlinear localized 
modes and consider the case of a power-law nonlinearity in detail. We discuss several scenarios 
of the instability-induced dynamics of the nonlinear impurity modes, including the mode decay or 
switching to a new stable state, and collapse at the impurity site. 



I. INTRODUCTION 



Wave scattering by localized impurities (or defects) is 
a fundamental problem of solid state physics [|l[. Im- 
purities break the translational symmetry of a physi- 
cal system and lead to several effects such as wave re- 
flection, resonant scattering, and excitation of impurity 
modes — spatially localized oscillatory states at the im- 
purity sites B. These two kinds of problems, i.e. wave 
scattering in inhomogeneous media and defect-supported 
localized modes, appear in many different physical prob- 
lems, such as the scattering of surface acoustic waves 
by surface defects or interfaces [Q, excitation of defect 
modes in superconductors in the vicinity of the twin- 
ning planes Q and high- Tc superconductors |^], the dy- 
namics of the tight-binding Holstein-type models of the 
electron-phonon coupling [pln|, light propagation in di- 
electric super-lattices with embedded defect layers [pj, 
excitation of defect states in photonic crystal waveg- 
uides Pi, light trapping and switching in non- uniform 
waveguide arrays [ [lO[ , etc. In all such cases, the impu- 
rities (or defects) lead to the energy trapping and local- 
ization in the vicinity of the defects, that occurs in the 
form of spatially localized impurity modes. 

When nonlinearity becomes important, it may lead to 
self-trapping and energy localization even in a perfect (or 
homogeneous) system in the form of intrinsic localized 
modes. Spatially localized modes of nonlinear systems 
are usually associated with solitary waves (or solitons) in 
continuous models, or discrete breathers in lattice mod- 
els; and they have been a subject of intensive studies dur- 
ing the past years |11|. However, the study of nonlinear 
phenomena in inhomogeneous and disordered systems is 
still largely an open area of research |12| . Simultaneous 
presence of nonlinearity and disorder is associated with 
various dynamical processes in solids, biological systems, 
and optics |l3| . For example, nonlinearities can become 
important even in a harmonic lattice due to the interac- 
tion of an exciton with the lattice vibrations [Q , where 
impurities may appear as a result of doping of materials 
with atoms or molecules that have stronger local cou- 



pling. Such impurities can also appear in (inherently 
nonlinear) spin- wave systems due to, e.g., a local varia- 
tion of the coupling between neighboring spins pq ] . 

When both nonlinearity and disorder are present si- 
multaneously, it is expected that competition between 
two different mechanisms of energy localization (i.e. one, 
due to the self-action of nonlinearity, and the other one, 
due to localization induced by disorder) will lead to a 
complicated and somewhat nontrivial physical picture of 
localized states and their stability. In this paper, we con- 
sider one of the examples of such a competition, and ana- 
lyze different types of nonlinear localized modes and their 
stability in the framework of the generalized nonlinear 
Schrodinger (NLS) equation with a point-like impurity. 

The problem we analyze here has a number of impor- 
tant physical applications ranging from the nonlinear dy- 
namics of solids [0 6| Jl4|Jl5t] to the theory of nonlinear pho- 
tonic crystals [ppj and waveguide arrays |10| in optics. 
In application to the theory of electromagnetic waves, 
this problem describes a special case of a stratified (or 
layered) dielectric medium for which nonlinear guided 
waves and their stability has been analyzed during the 
last 20 years ]l6| , p7t . For other apphcations, the theory 
of nonlinear localized modes is less developed and, in par- 
ticular, only a few publications |7||lq-p(]t addressed the 
important issue of stability of nonlinear impurity modes. 
In this paper, we study the properties of nonlinear im- 
purity modes in the framework of the generalized NLS 
equation and develop, for the first time to our knowledge, 
a systematic classification and linear stability analysis of 
spatially localized impurity modes of three distinct types. 
Our results can be linked to different special cases of the 
theory of nonlinear guided waves in layered dielectric me- 
dia, and they also provide a generalization of the theory 
of nonlinear impurity modes in solids, together with the 
analysis of their stability and instability-induced dynam- 
ics, emphasizing the cases where we can observe a clear 
evidence of competition between the two physical mech- 
anisms of energy localization. 

The paper is organized as follows. In Sec. we discuss 
our physical model and describe three different types of 



nonlinear localized modes supported by a point-like im- 
purity: symmetric one - an d two-hump modes and asym- 
metric modes. Section III includes a summary of the re- 



sults of the linear stability analysis, and it also presents 
the analytical criteria of the mode stability, for a general 
form of the NLS equation with a nonlinear impurity. A 
detailed analysis of the mode structure and stability, as 
well as the discussion of the competition between two dif- 
ferent mechanisms of energy localization, are presented 
in Sees. IV -VI, for the particular case of the power-law 



nonlinearity and both attractive and repulsive impuri- 
ties. At last. Sec. VII discusses two different types of the 
nonlinearity-induced collapse dynamics of the nonlinear 
modes, including collapse at the impurity site. 



nonlinear medium with a Kerr or non-Kerr response (see, 
e.g., Ref. |16|); in this case the time variable t stands for 
the propagation coordinate along the layer, and x is the 
transverse coordinate. 

For spatially localized solutions, Eqs. (|l|)-(g) conserve 
the power. 



P 



and the Hamiltonian 
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T{I';x)dI' \dx. 



(3) 



(4) 



II. MODEL AND LOCALIZED MODES 

We consider a general problem in which the dynam- 
ics of elementary excitations of a physical system (e.g., 
phonon, magnons, etc.) is described by an effective equa- 
tion for the wave-packet envelope ^{x^t) |21|. When 
the density of such quasi-particles becomes high enough, 
their interaction should be taken into account, e.g. in 
the framework of the mean-field approximation. In the 
simplest case, the quasi-particle interaction and collec- 
tive phenomena in an inhomogeneous medium can be 
described by the nonlinear Schrodinger (NLS) equation 
for the wave-packet envelope i}}{x,t)^ 



,dtp d'^tp 



+ 7^ + -^(^;^)^ = 0' 



dt dx^ 



(1) 



where / = jVP characterizes the density of the quasi- 
particles, t is time, x is the spatial coordinate, and the 
real function J^{I; x) describes both nonlinear and dis- 
ordered properties of the medium. We consider the case 
when the inhomogeneity is localized in a small region. 
Then, if the corresponding wavelength is much larger 
than the defect size, in the continuum limit approxi- 
mation the inhomogeneity can be modeled by a delta- 
function and, therefore, we can write 



T{I;x)^F{I)+6{x)G{I), 



(2) 



where the functions F{I) and G{I) characterise the prop- 
erties of the bulk medium and impurity, respectively. 
Hereafter we assume that the nonlinear term F{I) does 
not include a constant linear part, i.e. -F'(O) = 0, as oth- 
erwise it is always possible to re-scale the original Eq. (|l|) 
by introducing a new function ip = ipexp[—iF{0)t\. 

The model (|l|),(H) appears in different physical prob- 
lems of the macroscopic nonlinear dynamics of solids and 
nonlinear optical systems. In particular, it describes a 
special case of a more general problem of the existence 
and stability of nonlinear guided waves in a layered (or 
stratified) dielectric medium, where the delta- function 
defect corresponds to a very thin layer embedded into a 



Note, however, that the total momentum is not conserved 
since the translational invariance of the model (0)-(0) is 
broken by the presence of the inhomogeneity. 

We look for spatially localized stationary solutions of 
Eqs. (|l|)-(y) in the standard form, 

V(a;,t) ^u{x)e"^\ 

where w > is the mode frequency (w > for the mode 
to be exponentially localized), and the real function u{x) 
satisfies the equation: 



d^u 
dx'^ 



+ -r^+ F{I)u + 5{x)G{I)u = 0. 



(5) 



Let us first discuss the well-known case when the im- 
purity is absent, G{I) = 0. Then, a localized solution 
uo{x) of the reduced Eq. (||), 



LOUf) + 



Uq 



dx-^ 



F{ul)u^ = Q, 



(6) 



describes a self-trapped state in a uniform nonlinear 
medium. Due to the translational symmetry, a local- 
ized solution of Eq. (g) can be presented in the form 
^0(3; — 3^0)1 where xq is an arbitrary position. We also 
note that the mode profile is symmetric, uq{x) = uq{—x), 
it does not contain zeros, uo{x) > 0, and it has a single 
hump since duo/d\x\ < at a; 7^ 0. Such a localized 
state is possible in a self-focusing medium in the form of 
a bright solitary wave (or soliton) , which is characterized 
by the power 



Po(^) 



Uq{x) dx. 



(7) 



In contrast, in a defocusing medium, spatially localized 
solutions of Eq. (|6|) do not exist. However, certain singu- 
lar solutions with zero asymptotic at a; ^ ±00 are still 
possible in this case, and they will play an important role 
for constructing localized impurity modes. 

When a point- like defect is introduced into the system, 
the translational invariance of the model is broken at the 
defect location, x = 0. Nevertheless, the nonlinear modes 



of the inhomogeneous model (g) localized near the impu- 
rity site X = can be constructed by using the solutions 
of the homogeneous equation (g). Indeed, such a solution 
should satisfy Eq. (^ on both sides of the defect, and it 
can therefore be presented in the following general form. 



u{x) = 



Uq{x - Xq), X > 0, 
Uq{x — SXq), X < 0, 



(8) 



where xq and s are yet unknown parameters. Our task 
now is to satisfy Eq. (g) at x = in order to define the 
unknown parameters xo and s through the impurity char- 
acteristics. The first constraint follows from the field con- 
tinuity at the impurity site, u{0~^) — u{0~). For Eq. (||) 
it immediately yields the condition s — ±1. Thus, the 
parameter s defines the symmetry of the localized mode, 
i.e. the mode is symmetric for s — —1 and asym,m,et- 
ric for s = -t-1. In order to derive the second matching 
condition, we integrate Eq. (p) over an infinitely small 
segment around the impurity point a: = and obtain the 
transcendental equation for the parameter xq: 



(1 - s) — {xo) = G{Io)uo{xo), 
- ulixo). 



(9) 



where Iq 

Let us discuss the general properties of the impurity 
modes as localized solutions of a homogeneous medium 
satisfying the matching conditions obtained above. Sym- 
metric modes (s = —1) can be presented in the form 
u{x) — uq{\x\— Xq), and have the power 



P{uj) = 2 



+ 00 



Uq(x) dx. 



(10) 



Using the properties of the homogeneous solitary wave 
solution uq{x) we see that if the impurity is attractive 
(i.e. G{Io) > 0) then xo < 0. Thus the resulting sym- 
metric profile u{x) is single-humped [see Fig. [l|(a)] and 
P < Pq. When the impurity is repulsive (i.e. G(/o) < 0) 
then Xq > and thus the resulting profile is double- 
humped with P > Pq, eiS shown in Fig. |l|(b). We note 
that the single-hump solution may still exist in a self- 
defocusing medium being constructed from two pieces of 
a singular solution of the homogeneous model (see be- 
low), whereas the two- hump solutions are possible only 
in a self-focusing medium. 

For asymmetric modes s = -t-1 and the left-hand side 
of Eq. (||) vanishes. Thus solutions are possible only if 
the impurity response vanishes as well, i.e. G{Io) — 0. 
This condition can be satisfied when the nonlinear and 
linear parts of the impurity response compensate each 
other, i.e. for certain values of Iq. Remarkably, the pro- 
file of the asymmetric impurity mode coincides with that 
of a bright soliton in an uniform self-focusing medium, 
u{x) = uo{x — Xo), and thus P = Pq. However, its po- 
sition is fixed by the impurity according to the relation 
G(/o) = [see Fig. ||(c)]. Naturally there always ex- 
ist two degenerate asymmetric solutions corresponding 
to positive and negative xq, respectively. 




FIG. 1. Characteristic profiles of nonlinear impurity 

modes in a self-focusing medium: (a,b) symmetric modes sup- 
ported by an attractive {xo < 0) and repulsive (xo > 0) impu- 
rity, respectively; (c) asymmetric impurity mode with xq > 0. 



III. STABILITY ANALYSIS 



A. General formalism 



One of the key properties of a nonlinear localized mode 
is its linear stability determined by the character of the 
mode dynamics under the action of small perturbations 
of its stationary state. In general, two different scenarios 
of the perturbation-induced mode dynamics are possible. 
In the first case, a nonlinear mode can acquire only small 
distortions to its steady-state profile and the parameters 
of a nonlinear mode oscillate in the vicinity of its sta- 
tionary state. In this case, we call such a nonlinear mode 
linearly stable. On the other hand, under the influence 
of small perturbations initial deviations of the nonlinear 
mode parameters from their stationary values can grow 
exponentially; and in this case we define the nonlinear 
mode as linearly unstable. 

We consider localized modes that are square- or L^- 
integrable (have finite power or L^-norm) , and have pro- 
file functions that belong to a Hilbert space whose inner 
product is the L^-norm 



{a{x),b{x)} 



a*{x)b{x) dx, 



(11) 



where a*{x) denotes the complex conjugate of the func- 
tion a{x). To find the hnear stabihty conditions we 
consider the evolution of a small-amplitude perturbation 
f{x,t) of the stationary solution, i.e. 

ijix,t)^[uix) + fix,t)]e'^K (12) 

Writing the perturbation in the form 

fix,t) = [v{x) - wix)] e'"* + [v*ix) + w*{x)] e-'^'\ 

(13) 

and substituting Eq. ( p^ into Eq. (P, we obtain an 
eigenvalue problem for the functions vyx) and 'w{x), 



Lqw — Qv, Liv — Qw, 
where f2 is the complex spectral eigenvalue, and 

L^ = -(f/dx^ +Lu-Uj, Uj = Fj + Gj5{x) 
Fo = F{u^), Fi = Fo + 2u^F'iu^), 
Go = G(/o), Gi = Go + 2/oG (/q). 



Here prime denotes differentiation with respect to the ar- 
gument. It follows from Eq. (|l3|) that the mode is stable 
if all the eigenvalues fi are real and unstable otherwise. 
To proceed we reduce Eqs. (nih to a single equation: 



LqLiv = i^ V, 



(15) 



where stability requires all eigenvalues fl^ to be positive. 
It is straightforward to show that Lq = L^L^, where 
L = zLd/dx + u~^{du/dx), and thus one can instead 
consider the auxiliary eigenvalue problem, 



L^LiL+i 



= n^v, 



which reduces to Eq. ( p^ ) after the substitution v = L+-C. 
Because the operator i+ has only a single neutral mode, 
v{x) = u^^{x), which is not an eigenmode of Eqs. ( [15[ - 
Uw , these two eigenvalue problems have equivalent spec- 
tra (see Ref. |20)). Since the operator L~ LiL^ is Hermi- 
tian all eigenvalues Vl^ of Eqs. dlq) and (16) are real. 

The operators Lj {j = 1, 2) are well-studied in the lit- 
erature, in particular as a characteristic example of the 
spectral theory of second-order differential operators (see 
Ref. [p2[). For our problem we use two general mathemat- 
ical results about the spectrum of the linear eigenvalue 

(i) the eigenvalues can be or- 



problem Ljip]i' = \n\>n^ 



dered as A, 



(i) 

Tl + l 



> 



x(i) 

An 



where n > defines the number 

U). 



of zeros in the corresponding eigenfunction (pn \ (ii) for 



Let us fist discuss the properties of the operator Lq 
whose spectrum we define as A„ . The neutral mode of 
Lq is the stationary solution of Eq. (||), since Lqu{x) = 0. 
Moreover, u{x) > is the ground state solution with no 



for n > 0. This 



nodes and therefore A„ > Aq 

means that the operator Lq is positive definite on the 
subspace of functions orthogonal to u{x), which allows to 
use several general theorems p^ , p3| -p6[ to link the sta- 
bility properties to the number of negative eigenvalues 
of the operator Li, whose spectrum we define as A„ : 
(i) the mode is unstable if there are two (or more) neg- 
ative eigenvalues, i.e. if Aj^ < 0; (ii) the mode is stable 
if Li is positively definite, i.e. if Aq > 0; (iii) in the in- 
termediate case the stability depends on the slope of the 
function Pj^) according to the Vakhitov-Kolokolov (VK) 
criterion [2^, i.e. the mode is stable if dP/du > and 
unstable otherwise. Thus, to distinguish between these 

(14) cases, it is suflficient to determine the signs of Aq and 

In general, the spectral properties of the linear opera- 
tor Li depend on the mode frequency u, and the number 
of its negative eigenvalues can change. This is associated 
with the so-called critical points in the power dependence 
P{uj) for the localized mode |]1^. One type of such crit- 
ical points is a bifurcation into two families of solutions, 
i.e. symmetric and asymmetric ones, which exist for the 
same mode frequency. Thus, the stability properties of 
these distinct types of nonlinear localized modes can be 
expected to be different, so we discuss them separately. 



B. Symmetric modes 

For symmetric modes with s = — 1 the operator Li has 
the symmetric neutral mode with zero eigenvalue 

Li{du/d\x\) = 0, 



(16) for the special value of the defect response 



Gl' ^ 2 



^ duo{x 
. dx 



-1 ^2 



d'^UQ{x) 
dx'^ 



(17) 



For xo > the profile u{x) has two humps and the neutral 
mode corresponds therefore to the second eigenmode (^3 
(two nodes) with A2 — 0. Thus the two-hump symmet- 
ric modes are unstable (two negative eigenvalues). For 
xo < the neutral mode is the nodeless ground-state 



V6 



(1) 



with A, 



(1) 



0. Thus the single-hump symmetric 



modes are stable (no negative eigenvalues), if Gi = G" . 
For other values of the defect response Gi it follows 
from the spectral theorem that the eigenvalue of the neu- 
tral mode decreases when Gi > G" (deeper well) and 
increases when Gi < G". Therefore, for xq < and 



"deeper" potential wells, Uj{x) > Uj{x), the correspond- Qi < Gf we have a[,^^ > and the modes are stable. On 



ing set of eigenvalues is shifted "down', A„ < A 



A3) 



the other hand, due to the symmetry of the potentials in 



the linear eigenvalue problem, Uj{x) = Uj{—x), the am- 
plitude of the first-order eigenmode vanishes at the defect 
location (x = 0) , and thus its eigenvalue A^ does not de- 
pend on Gi. However, the inequality Aq < Aj^ < Aj 
is always fulfilled for any Gi, and thus A^ ^ < if xq > 
and Af ^ > if xo < 0. 

In the special case when xq = it is straightforward 
to show that the function du{x)/dx, which has a single 

zero at a; = 0, is a neutral mode of Li with Aj^ = 0. In 
Table | we summarize the properties of the operator Li 
and the corresponding general conditions for the stability 
of nonlinear localized modes. 

TABLE I. Stability conditions for symmetric modes 



C. Asymmetric modes 



The profiles of asymmetric nonlinear modes coincide 
with those of solitons (since Go — 0), but the spectrum 
of the operator Li can become different. Only in the 
special case Gi = 0, there exists a first-order (one node) 
neutral mode with Aj^ = 0, i.e. Li{du/dx) = 0. Then, 
from the spectral theorems [Q it follows that A]^ < 

if Gi > and A^^^ > if Gi < 0. Thus the asymmet- 
ric mode is always unstable for Gi > 0, with respect to 
translational shifts along the x axis. 



(1) 



Condition 




I/i spectrum 


Stability 


xo > 




Af,^' < X['^ < 


unstable 


xo < 0, Gi 


<Gf 


o<x\^Ux['> 


stable 


xo < 0, Gi 


> Gf , dP/duj > 


\^^ < < X['^ 


stable 


xo < 0, Gi 


> Gf , dP/duj < 


a^^^<o<aw 


unstable 



To determine the sign of Aq for Gi < we consider, 
without a lack of generality, the case xq > and study 
the operator Li with the potential 



Uiix) 



Fi [ul{x - xq)] , X > xq, 
— oo, X < Xo. 



(18) 



It is important to connect the spectral characteristics 
of the operator Li, and thus the stability properties, 
with the character of the power functional P{(jj)- To 
do so we notice that at the bifurcation pointy defined as 
P{lo) — Po{uj), the parameter xo and thus the eigenvalue 

Xi changes sign. The two-hump symmetric modes with 

P(w) > Po(w), for which xq > and thus A^ ' < 0, are 
therefore always unstable. Instead a new family of asym- 
metric modes emerges at the bifurcation point, whose 
stability properties we discuss in the next section. 

The single-hump symmetric modes (i.e. xo < 0) can 
only change their stability properties through a change 
of sign of the lowest eigenvalue Aq or the slope dP/diiJ, 



It is straightforward to check that the lowest eigenvalue 
of Li is zero with the ground state 



4^^ 



duo{x — Xo)/dx, X > Xq, 
0, X < Xq. 



For Gi < we have Ui (x) > Ui{x), and thus according to 
the spectral theorem the eigenvalues of the correspond- 
ing operators are related as A„ < A„ . Therefore, Ag is 
always negative, meaning that the VK theorem applies. 
We co mbine these findings with the general stability cri- 
and summarize the stability conditions 



teria |7| 
in Table H 




J- i, rn 1 1 Fl T J i- c J i-1 -i- 1 -4- TABLE II. Stability conditions for asymmetric modes 

accordmg to lable |1[ In order to hnd the critical points -^ •' 



associated with the transition of Aq through zero we 
employ the approach developed in Ref. |1^ and differen- 
tiate Eq. (H) with respect to the mode power P, treating 
cj as a function of P. This gives us the relation 

u{x){duj/dP) = -Li{du/dP). 

The points where dio/dP — or u = are thus critical 
points, since there the operator Li has a zero eigenvalue 
Ag — with eigenfunction du/dP. It is possible to show 
that these are the only critical points |0|. Moreover, it 
can be demonstrated that Gi > G" for the branch orig- 
inating from the critical point with the positive slope, 
dP/dijj > 0, while Gi < Gf if dP/duj < in the vicin- 
ity of the critical point. Therefore, the localized modes 
are always stable close to such critical points. This con- 
clusion implies, in particular, stability of modes with a 
vanishing power, P — > 0. In this linear limit the impurity 
can support a localized mode only when it is attractive 
(i.e. Go > 0), and the mode frequency is ujo — Go/4:. 



Condition 


Li spectrum 


Stability 


Gi >0 


a:,^'<a(^'<o 


unstable 


Gi < 0, dPo/dio > 


A^^'<0<Ai^' 


stable 


Gi < 0, dPo/dio < 


X^'> < < Af ' 


unstable 



We note that, for a given defect response G(/), there 
can exist several families of asymmetric modes each char- 
acterized by the intensity at the defect Iq, satisfying the 
relation G(/o) = 0. There are always at least two solu- 
tions with the same /q but ±xo. All these degenerate 
families have the same power and are thus indistinguish- 
able on the P{uj) diagram. To determine the stability 
we look at the bifurcation from the symmetric modes at 
P{o-!) = Po{oj) (or Xo = 0). Performing the analysis simi- 
lar to that in Ref. ||l^ , we find that after the bifurcation 
Gi < if the branch for symmetric modes is above that 
of the asymmetric modes, P{uj) > Po(w). In the oppo- 
site case Gi > and the asymmetric modes are always 
unstable. 



IV. POWER-LAW NONLINEARITY 

We now demonstrate the characteristic existence and 
stabiHty features of nonhnear locahzcd impurity modes, 
by applying the results of Sees. Q and III to the illustra- 
tive case of power-law nonlinearities. In this Section we 
consider the analytically obtainable results. The detailed 
numerical results are presented in Sees. M and VI. 



A. Solitons in homogeneous media 

To construct the general solutions we consider first the 
special case of a homogeneous medium with G = and 
the bulk power-law nonlinearity 



F{I) = pi", CT > 0. 



(19) 



For self- focusing bulk media p = +1 and for defocusing 
bulk media p = — 1 . The profiles of the spatially localized 
stationary solutions can be found in the form 



uo{x) = Aq 



cosh "'"' '^ {(jy^x) 
sinh^ '"^ {(Ty/ujx) 



P 
P- 



-1, 



-1, 



(20) 



i1/2ct 



where Aq = [(1 + ij)uj] ' . Solution ( |20| ) for p = +1 
is the well-known bright soliton of the generalized NLS 
equation that exists for any positive value of a. The so- 
lution for p = — 1 is singular, since bright solitons do not 
exist in a self-defocusing medium. This singular solution 
will be used below to construct the profiles of the (regu- 
lar) localized modes in the general case with an impurity. 
To find the soliton power in a self-focusing medium 
(p — +1) we substitute the corresponding solution (|2C 
into Eq. (0) and obtain the result: 



Po{uj) = LU 2" 



22/^(1 -ffT)V'^r|(i/g) 

2a r5(2/fT)^ 



(21) 



where Ts is the standard Gamma-function. Applying the 
VK theorem |23| we find that the solitons are stable for 
a < O'er = 2, because in this case dPo/diu is positive. For 
(7 > 2 the solitons are unstable {dPo/dtu < 0) and they 
either collapse or decay P3'E3 ■ 



B. Nonlinear impurity in a linear medium 

We now study the effects due to defect-induced local- 
ization only, i.e. we consider the special case when a 
nonlinear defect with the power law-response 



G{l) = a + I3r, 7>0, 



(22) 



is embedded into a linear medium with F = 0. Chang- 
ing the signs of the parameters a and /3 we can describe 
both linearly and nonlinearly attractive and/or repulsive 



impurities. Assuming that a y^ this parameter can be 
rescaled to a = ±2, which we use in the following. 

As bright solitons do not exist in a linear medium, 
asymmetric modes are not possible as well. Using 
Eqs. (||),(||), and (||) we obtain the spatial profile of the 
symmetric one-hump localized modes 

u,{x) = A,eM-Mx\), Ap = {2V^~a)/P, (23) 



and the corresponding power. 






1/7 



(24) 



Let us analyze the existence properties of these modes: 
For a > a linear impurity mode appears at the cut-off 
frequency loq = a'^/A. If the nonlinearity is attractive 
{(3 > 0) then a whole family of localized modes exist 
with frequencies above cut-off (u > loq). If the nonlin- 
earity is repulsive {(3 < 0) then the family exists below 
cut-off (0 < w < ujq)- For a < no localized modes 
exist for repulsive nonlinearities (/3 < 0), whereas attrac- 
tive nonlinearities (/? > 0) supports localized modes at 
all frequencies a; > 0. 

To analyse the st ability properties we follow the ap- 
proach of Sec. IIIB and calculate the defect response 



G, = Gl'- + 2(3jA^\ 



(25) 



where the critical value G" = 2\/aJ is defined from 
Eq. (O). From the general results summarized in Ta- 
ble ffl, we see that the localized modes are always sta- 
ble when the nonlinearity is repulsive (/3 < 0), since 
Gi < G". When the nonlinearity is attractive (/? > 0), 
Gi > Gf and the VK criterion applies. From Eq. ( |24| ) 
we obtain that the sign of dPi/dco is given by the sign of 
[2(1 — 7)v^ + aj]. Thus we identify the critical power 
7cr = 1, which is half that in a homogeneous medium. 
For a > (and (3 > 0) the modes are therefore always 
stable for sub-critical powers 7 < 7cr. For 7 > 7cr high 
frequency modes with lu > uji = ll)q/{1 — 7cr/7)^ are 
unstable, whereas low frequency modes [uj < loi) are sta- 
ble. The opposite occurs for a < (and (3 > 0): In this 
case only high frequency modes {uj > uii) are stable for 
7 < 7cr, whereas low frequency modes are unstable. For 
powers above the critical value, 7 > 7cr, all modes are 
unstable. 



C. Nonlinear impurity modes 

We now consider a more general case when both the 
bulk medium and the defect have power-law nonlinear 
responses, as defined by Eqs. ( |l^ ) and (p2|). 

Symmetric nonlinear modes: Substituting Eq. (pfl) 
into Eq. (O) we obtain the relation: 



2pLU 



1/2 



f3iu^il 



af\l-pr=0, 



(26) 



which defines the spatial shift a;o for all values of the 
mode frequency lo. Here T = 7/cr and 



p{uj) = 



tanh (croj-'^'^xo), p 
coth ((Tw^/^xo), p 



-1, 
-1. 



(27) 



We see that for \p\ < 1, the solutions of Eq. ( [26| ) corre- 
spond to the modes in a self-focusing medium (p = +!)• 
Note that localization in a defocusing medium {p = —1) 
is only supported by an attractive impurity (Gq > 0), so 
that xo < and p < —1. In both cases the linear limit 
xq — > — 00 corresponds to p —f —1. Impurity response for 
a localized mode can be expressed in terms of the new 
variable, 

Go — -~2py/uj, 

Gi = Go + 2/?7/o^ = (27 + l)Go - 27a, (28) 

GT^-2V:j[p{l + a)-a/p], 

which will be useful in the following analysis. 

In order to determine the stability of the localized 
modes, we have to construct the power vs. frequency 
diagram P{uj). After a change of variables in Eq. ([l(]| ) we 
obtain the general expression for the power 



P{uj) = UJ 



■!-<, 2fl+g) 



1/<T 



|i-y 



2|-l + l/cr 



dy 



(29) 



To separate the effects induced by the bulk and the im- 
purity we re-write this expression in the form: 



P(co) = Po(^)^(p;^), 



(30) 



where the soliton power in homogeneous media, P^iiS)^ is 
defined by Eq. (|2|). Although bright solitons do not exist 
in defocusing homogeneous media, Eqs. ( |2l|) and ( pO|) can 
still be used for both types of bulk nonlinearity, i.e. for 
p = ±1. The functional ^(p, u) changes monotonically 
with the parameter p, so that p{d^/dp) > 0. In the linear 
limit when a > and uj —> ujq the mode power vanishes 
and thus ^(— 1;ct) = 0. Furthermore, for a self-focusing 
bulk medium (p = +1) the identity ^(— p; a) = 2 — ^(p; a) 
holds. Thus C(0;cr) = 1 and ^{l;a) = 2. The functional 
^{p; ct) cannot be expressed in elementary functions, but 
must be calculated numerically. 

Let us outline some general properties of the symmet- 
ric modes: Consider first linearly attractive impurities 
(a > 0) and the low-intensity limit when the mode fre- 
quency is close to the cut-off loq. From Eq. (|2^) we find 
the nonlinearity-induced frequency shift 



UJ 
UJQ 



1 



/32r(l-H 
2(1 +rt 



aV Lj, 



^1/2 





r> 1. 



i+pf, r<i. 



(31) 



It follows that for r > 1 the sign of the shift depends 
only on the bulk nonlinearity: the frequency is shifted 
up in self- focusing bulk media {p = +1) and down in 



self-defocusing bulk media (p = —1). For F < 1 the shift 
depends only on the defect nonlinearity: the frequency 
is shifted up when /3 > and down when (3 < 0. How- 
ever, Eq. ( |3l| ) presents only an asymptotic result valid at 
vanishingly small intensities. A simple graphical analysis 
of Eq. (Pq ) demonstrates that competition of the defect 
and bulk nonlinearities for pf3 < can lead to a more 
complicated (e.g. multi-valued) structure of the power 
dependence P{uj). 

Secondly, we analyze the properties of high-frequency 
modes in self- focusing media {p = -1-1). Analysis of 
Eq. (^) reveals that such modes always exist. Substi- 
tuting approximate solutions for p(w) when w — > c» into 
Eq. (E9) we find the power of the high-frequency modes 

r Po(t^)(l-/9Cc.^^^/'), r<i/2, 

P(w) ^ I Pi{uj), F > 1/2 and /3 > 0, (32) 

[ 2Po(w), r > 1/2 and /3 < 0, 

where G is a frequency-independent positive coefRcient. 
It follows from Eq. (^2|) that for a nonlinearly repulsive 
defect (/3 < 0) we have P{uj) > Po(w) and thus the high- 
frequency symmetric modes have two humps and are un- 
stable. For nonlinearly attractive defects (/3 > 0) the 
modes have one hump, i.e. — 1 < p < 0. It is straightfor- 
ward to check using Eq. (gsj) that under these conditions 
Gi > Go > Gf , and therefore according to Table | the 
mode stability follows directly from the power slope. In 
particular, if F < 1/2 then the nonlinearity of the bulk is 
effectively stronger than that of the defect, so that the lo- 
calized modes resemble the solitons in homogeneous me- 
dia, which are stable for a < 2. In the opposite case, 
F > 1/2, the modes resemble the nonlinear defect modes 
of linear media, which are stable for 7 < 1. Combin- 
ing these results we come to the conclusions: Symmetric 
localized modes supported by a nonlinearly attractive im- 
purity (/3 > 0) in a self-focusing bulk medium are stable 
in the high-frequency limit if <t < 2 and 7 < 1 simulta- 
neously. 

Finally, we study the properties of high-frequency lo- 
calized modes in self-defocusing bulk media {p = —1). 
In this case Eq. (|2^) has solutions for w — > -Fcx) only for 
a "strongly" nonlinear attractive impurity with F > 1/2 
and /? > 0. The corresponding modes have the same 
properties as in a self-focusing bulk medium, because 
the bulk nonlinearity acts only as a small perturba- 
tion. Indeed, according to Eq. (28) we have Gi — Gf — 
—2y/uJ[a/p-\-p{2j—a)] — 2-fa > for uj ^ 1, and therefore 
the mode stability follows from the power slope. 

Asymmetric nonlinear modes : Asymmetric modes can 
exist in self-focusing bulk media when the impurity re- 
sponse vanishes, G(/o) = 0. This is possible when a/3 < 
and the mode intensity at the defect site is /q = |a//?|^/^. 
Then, the spatial shift xq is given by the relation 



cosh(crVcJa;o) — A^/Il 



r/2 



(33) 



Asymmetric modes bifurcate from the symmetric modes 
at the bifurcation point Wb = /q /(I -h cr) when Aq = Iq 



and they exist for higher frequencies only. An asymmet- 
ric mode is actually a soliton trapped by the defect and 
is therefore only stable for cr < 2. However, the trap- 
ping (defined by the condi tion G gjlo) — 0) is stable if 
Gi = 2/37/0 < (see Sec. IIIC ), which is satisfied for 
nonlinearly repulsive impurities with f3 < 0. 



V. SELF-FOCUSING MEDIUM 



In this section we present a detailed numerical analy- 
sis of the nonlinear impurity modes for power-law non- 
linearities in a self-focusing bulk medium {p — +1). The 
self-defocusing bulk medium is treated in Sec. VI. We 



term the impurity attractive when it supports a localized 
mode in the linear limit (a > 0) and has a self-focusing 
(attractive) nonlin ear r esponse (/3 > 0). This case is 



considered in Sec. V A . The same impurity, but with 
a self-defocusing nonlinear response {j3 < 0) is ter med a 
m,ixed impurity that supports linear modes (see Sec. |VB| ). 
The case when the defect is repulsive in the linear limit 
(a < 0), but may become attractive for larger intensi- 
ties (/3 > 0) due to its nonlinear properties is the most 
complicated. This mixed impurity with attractive nonlin- 



earity is investigated in Sec. V C, When both a and (3 are 



negative the properties of localized modes in self-focusing 
bulk media are trivial: they always have two-hump sym- 
metric profiles and are therefore always unstable. 




FIG. 2. Power vs. frequency diagram for 7/(7 — 1 and 
two characteristic mode profiles corresponding to the points 
(a) and (b). Dotted curve: power Pi{u}) of impurity modes in 
linear bulk media. Parameters: a — 1, ^ = 1, a — 2, l3 — 1. 



We first consider the influence of defect nonlinearities 
that enhance the attraction (/3 > 0). In this case the 
asymmetric modes do not exist. For the symmetric lo- 
calized states we always have xq < and, as follows from 
Eq. (11), Gi > Go > GY, so that the VK criterion |2|] 
applies according t o Tab le ||. To illustrate the general 
results given in Sec. IV C, we present examples in Figs. 
and H, and the stability diagram in Fig. ^. Localized 
modes always exist above cut-off, uj > luq, and the func- 
tional P{ui) is single- valued. Such a behavior is observed 
because both the defect and the bulk medium have self- 
focusing nonlinearities, which induce a positive frequency 
shift, as mentioned earlier. 



A. Attractive impurity 



Let us consider self- focusing bulk media {p — +1) and 
defects that support localized modes in the linear limit 
(a > 0). In this case the symmetric mode branch on the 
P(w) diagram always starts at w = loq^ where P vanishes. 
As follows from the stability analysis this point is critical, 
and therefore the localized modes are always stable in its 
close vicinity, see Figs. and pi However, for higher fre- 
quencies Lo the properties of the localized modes depend 
on the characteristics of the nonlinearity, determined by 
the powers 7 and a and the coefficient /3. 



The parameters for Figs. ^ and ^correspond to quali- 
tatively different cases. For 7/cr > 1/2 (Fig. 0) the defect 
nonlinearity is effectively stronger than the bulk nonlin- 
earity, and thus the modes resemble the nonlinear defect 
modes in linear bulk media. Since 7 and a do not ex- 
ceed their critical values 7cr = 1 and Ccr = 2 [see Fig. |^, 
point (a)], all modes should be stable. For -f/a < 1/2, 
the high-frequency symmetric modes resemble the soli- 
tons of homogeneous media [see Fig. §(b)]. Therefore, 
these modes are unstable if the bulk nonlinearity is super- 
critical a > (Tcr [see Fig. H. 
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FIG. 3. Power vs. frequency diagram for y/a = 1/4. Solid 
and dashed thick curves show stable and unstable modes, re- 
spectively. Thin dashed curve: soliton power Po{u}). Thin 
dotted curve: power Pi{iL}) of impurity modes in linear bulk 
media. Modes corresponding to points (a) and (b) are shown 
below. Parameters: a = 4, 7 = 1, a = 2, /? = 0.1. 
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FIG. 4. Stability regions for symmetric localized modes 
with Q > and /3 > in self-focusing media {p — -1-1): I: the 
modes are stable and exist for lj > o^o; II: stable modes exist 
only near cut-off cj ~ ljq- Points (a) and (b) correspond to 
the cases shown in Figs. S and M, respectively. 



To study the instability-induced dynamics of symmet- 
ric localized modes we solve Eq. (M numerically with 
slightly perturbed modes as initial condition. Depending 
on the perturbation, two different scenarios are possi- 
ble. If the power is initially decreased, then the mode 
spreads out and transforms into a lower-frequency stable 
mode, as shown in Fig. |g(a). We note that this switching 
process is accompanied by some power loss due to radi- 
ation. An initial increase of the mode power can lead 
to collapse if the nonlinear self-focusing dominates the 
linear diffraction. Figure ^(b) shows a collapsing mode, 
whose amplitude goes to infinity in a finite time due to 
the effect of a supercritical bulk nonlinearity, cr > 2, even 
though the collapse occurs at the impurity site. This 
collapse instability was earlier investigated for uniform 



nonlinear media [g^ , with many of the general features 
applying to the impurity modes as well. For example, the 
collapsing solution typically consists of a slowly evolving 
background and a highly-localized central part having an 
almost self-similar profile. 




(a) 




(b) 




FIG. 5. Evolution of an unstable mode [Fig. H(b)] with 
the input power decreased (a) or increased (b) by 0.1%. Two 
different scenarios of the instability-induced mode dynamics 
are observed: (a) switching from an unstable to a stable state, 
and (b) collapse induced by the bulk nonlinearity. 



B. Mixed impurity that supports linear modes 

We now consider defects that have a repulsive nonlin- 
earity, /3 < (and still a > 0). Due to the competition 
between focusing bulk and defocusing defect nonlincar- 
ities the power functional of the symmetric modes can 
become multi-valued, with two or three states having the 
same frequency (see Figs. ^ and |^). Such states can even 
exist below the linear cut-off frequency, i.e. for w < wq. 
This is always possible for 7 < a, because the branch 
starting at the critical point P = has a negative slope 
(still stable since Gi < Gf), and then the slope changes 
sign as the branch goes through another critical point, as 
shown in Fig. 0. For 7 > cr the initial slope is positive, 
but then the branch can go through two critical points, 
which appear if /3 < /3cr (see Fig. 0), where 



(3c 



,i-2r (r-1) 



(2r - 1) 



2(2r-i) 



(i + a)r2(r-i) 



(34) 




FIG. 6. Power vs. frequency diagram for 7 > cr (a = 1, 
7 = 2, a = 2, /3 = -1, /?cr ^ -0.21). Solid and dashed thick 
curves show stable and unstable localized modes, respectively. 
Open circle: bifurcation point from symmetric (S) to asym- 
metric (AS) modes. Thin dashed curve: soliton power Po('^). 
The modes corresponding to points (a), (b), and (c) are shown 
below. 
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FIG. 7. Power vs. frequency diagram for 7 < cr (ct = 2, 
7 = 1). Parameters and notations are the same as in Fig. |6| 
Dotted curve: impurity mode power in a linear medium, 

Rico). 

The branch of asymmetric modes emerges at the bi- 
furcation point where P{uj) coincides with the power of 
solitons in homogeneous media, Po{uj) (open circles in 
Figs, g and m. As demonstrated in Sec. IV C the stabil- 
ity of asymmetric modes for /3 < is the same as that 
of solitons in bulk media (see Sec. IV A), i.e. the modes 
are stable for a < 2 and unstable otherwise. After the 
bifurcation point, i.e. for uj > uJb, the symmetric localized 
mode becomes two-humped and therefore unstable. 

In Figs. g(a) and (b) we show the evolution of two- 
hump symmetric modes when the bulk nonlinearity is 
sub- and super-critical, respectively. In both cases a sym- 



metry breaking occurs, and the mode is repelled to one 
side of the defect. In case (a) a stable asymmetric state 
is excited (as u < 2), together with its internal mode, 
leading to a slowly decaying quasi-periodic beating. On 
the other hand, if bright solitons are unstable, the mode 
collapses in the bulk medium as shown in Fig. |8|(b). Note, 
that the amplitude at the impurity always remains finite 
due to the repulsive nonlinear response of the defect. 





FIG. 8. (a,b) Evolution of perturbed unstable two-hump 
symmetric localized modes shown in Figs. H(c) andM(c), re- 
spectively. In (a) the mode transforms into a stable asymmet- 
ric mode; in (b) the transition occurs to an unstable asymmet- 
ric mode, which subsequently collapses in the bulk medium. 



C. Mixed impurity with attractive nonlinearity 

We finally consider the most complicated case of a lin- 
early repulsive, but nonlinearly attractive defect (a < 0, 
P > 0). Some properties of the symmetric localized 
modes are revealed by the symmetry constraints follow- 
ing Eqs. (p^) and (pOh. Specifically, under an inversion of 
the defect response, {a,f3) -^ (— a,— /3), the normalized 
shift p{xa) and power P change as follows: p — > (1 — p) 
and P{lu) -^ 2Po{uj) — P{io). Applying this transforma- 
tion to the case a < and /3 > and using the results 
obtained above, we find that the branch of symmetric lo- 
calized modes in the Piuj) functional starts at the point 
P{u>o) = 2Pq. In the vicinity of this point the symmet- 
ric modes are unstable because xq > (see Figs. ^ |l^, 
and O). At the bifurcation point (open circle) the family 
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of asymmetric modes emerges, having the same power as 
sohtons in a bulk medium, Pq{uj). The symmetric modes 
have one-hump profiles after the bifurcation. 
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FIG. 9. Power vs. frequency diagram for a — 1, "/ = 0.5, 
P — 1, and a — —2. Notations are the same as in Fig. pi 
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FIG. 10. Power vs. frequency diagram for a = 1, 7 = 1.5, 
l3 = 0.01, and /?max — 0.0262. Parameters and notations are 
the same as in Fig. H. 



The stability of symmetric one-hump localized modes 
depends on the nonlinearity characteristics. Since we 
have /3 > and xq < 0, it follows from Eq. (^8|) that 
Gi > Gq > Gf\ and stability is determined by the power 
slope according to T able 0. In particular, as demon- 
strated in Sec. [VC, the modes are stable in a range 
of frequencies unlimited from above if cr < 2 and 7 < 1 
(see, e.g.. Fig. ^. This parameter region I in Fig. ^ is 
identical to that for a > (see Fig. 0); the mode prop- 
erties are also similar, because in the limit uj ^ luq the 
linear response of the defect acts as only a small pertur- 
bation. However, there is an important difference: while 
for a > the modes are always stable for u; ~ wo, this is 
not so for a linearly repulsive impurity. Although in the 



latter case the high-frequency modes are also unstable if 
the nonlinearity powers exceed critical values, we found 
that for the nonlinearity parameters corresponding to re- 
gions II or III in Fig. [ij the modes can still be stabilized 
at some frequencies. 
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FIG. 11. Power vs. frequency diagram for a — 2.5, 
7 = 0.5, /3 = 3.3, and /3min — 3.17. Parameters and nota- 
tions are the same as in Fig. H. 

The first mechanism of stabilization in a bounded fre- 
quency region above the bifurcation frequency tot is re- 
alized when the modes are close to stable solitons (for 
CT < 2) [see Fig. |l^, and the strength of the defect is 
relatively weak (region II in Fig. |2|). This is possible 
if the nonlinearity coefficient docs not exceed a thresh- 
old value, which can be calculated from the condition 



< ^ < ^max = 



(1 



(2-CT)22/'^r|(l/CT) 



47|a|r5(2/CT) 



2r 



(35) 



On the other hand, if ct > 2 the modes can be sta- 
bilized by the attractive defect nonlinearity, as demon- 
strated in Fig. O, provided the impurity supports non- 
collapsing highly localized states (i.e. for 7 < 1, see region 
III in Fig. n2), and has a sufficiently strong nonlinearity, 
/3 > /?min- The minimum value /^min cannot be deter- 
mined analytically. Its characteristic dependencies on ct 
for several values of 7 are shown in Fig. O. 

As demonstrated in Sec. IVCj the asymmetric modes 
are unstable for /3 > with respect to a translational 
shift along the x axis. To study the development of this 
instability we perform numerical simulations. The results 
confirmed that the mode is either attracted or repelled by 
the impurity, depending on the type of perturbation, as 
illustrated in Figs. |l^(a) and |l^(b), respectively. These 
examples correspond to the case of sub-critical defect 
and bulk nonlinearities, and thus the mode eventually 
transforms into a stable impurity mode or a moving soli- 
ton. However, for stronger nonlinearities the instability 
induced dynamics can result in a collapse. 
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FIG. 12. Stability regions of localized modes for impuri- 
ties with a > and /3 > in self-focusing bulk media: I: the 
modes are stable for frequencies above a threshold, lu > ujth', 
II and III: stable localized states can exist in a bounded fre- 
quency range. In the region without shading all modes are 
unstable. Parameter values at points (a), (b), and (c) corre- 
spond to Figs. H, hd, and hll respectively. 




FIG. 13. Parameter /3min versus bulk power nonlinearity 
CT for 7 = 0.1 (dotted), 0.5 (solid), and 0.8 (dashed). 




FIG. 14. (a,b) Evolution of the perturbed asymmetric 
mode shown in Fig. Mh). (a) Attraction by the impurity and 
transformation into a stable symmetric mode, (b) Repulsion 
by the impurity and transformation into a moving soliton. 



VI. SELF-DEFOCUSING MEDIUM 

In this section we focus on self-defocusing bulk media 
(p = — 1). In contrast to self-focusing bulk media investi- 
gated above, bright sohtons cannot exist in homogeneous 
self-defocusing bulk media, and thus localized modes ap- 
pear solely due to the presence of the impurity, having 
always a one-hump profile. Again we consider the differ- 
ent cases separately. 



A. Attractive impurity 



First we consider attractive defects with both a > 
and /? > 0. It follows from the analysis presented above 
that locahzed modes exist and are stable near the cut-off 
frequency, cj ~ wq- For defects with a purely linear re- 
sponse (7 -^ 0) there exists a single branch of localized 
solutions on the P{tu) diagram in the region < uj < luq. 
This case was considered in Ref. [M for a = 1. 

Remarkably, defects with a nonlinear response can 
support localized modes with frequencies above the lin- 
ear cut-of f frequency. Indeed, as we demonstrated in 
Sec. IV C, the initial slope of the P{uj) functional for 
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weakly nonlinear modes is determined by the nonlinear- 
ity of the impurity for F < 1. If F < 1/2 the branch 
goes through a critical point and disappears as w — > 
[see Fig. |l5|(c)], while the power P remains bounded for 
cr < 2, and is unbounded for cr > 2. For 1/2 < F < 1 
the critical point appears only for (3 < Pen where Per is 
defined by Eq. (M) [see Fig. [l5|(a)]. 




cussed above. The lower branch approaches the linear 
limit P{u'o) = and has initially a negative slope, which 
changes at the critical point for 1/2 < F < 1, as shown 
in Fig. |l^(a) (consistent with the inequality /3o < Pa-)- 




FIG. 16. Stability regions of localized modes for defects 
with Q > and /? > in defocusing media. I: modes ex- 
ist and are stable at both low and high frequencies; II: stable 
modes exists only close to and below the linear cut-off a; ~ loq. 
Points (a)-(c) correspond to Figs. Qq(a)-(c), respectively. 




FIG. 15. Power vs. frequency diagrams for (a) 

1/2 < r < 1 (a = 1, 7 = 0.7, P = 0.52, /?„ ^ 0.603, and 
Po ^ 0.533); (b) r > 1 (a = 1, 7 = 1.5, /3 = 0.11, and 
Po ~ 0.105); (c) r < 1/2 (a = 3, 7 = 1, /? = 0.5). Nota- 
tions are the same as in Fig. ti and a — 2. The marked point 
in (a) corresponds to an unstable solution whose instability 
dynamics is illustrated in Fig. hq. 

If the nonlinear response of the impurity is effectively 
stronger than that of the bulk, highly localized hig h fre- 
quency modes can exist. As was found in Sec. IV C 
this occurs for F > 1/2. These modes have the power 
P ~ -Pi(w) when lli ^ llIq, and are thus stable for 7 < 1 
(region I in Fig. |l^. The mode frequency is bounded 
from below, uj > lu nun, for P > Pq, where 



Po{a,a,r) = 



|2F-l|(l + a)i 



2F- 1 



aT 



2r 



(36) 



as shown in Fig. |l^(b). For P < Pq the functional P{u)) 
has two branches originating at u; = 0, with the upper 
branch corresponding to the highly localized modes dis- 



» 




FIG. 17. Evolution of a perturbed asymmetric impurity 
mode corresponding to the marked point in Fig. |l5|(a). (a) 
Switching to a higher-frequency mode of the same branch, (b) 
A reduction of the power by 5% leads to transformation to a 
broader low-frequency mode of a lower branch. 

To clarify some features of unstable modes in this case, 
we perform numerical simulations for the mode corre- 
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spending to the marked point in Fig. |l^(a). We find that 
all small perturbations result in switching to a more lo- 
calized state of the same (upper) branch, as illustrated in 
Fig. |l3(a). However, if the power is decreased below the 
minimum of the upper branch, the mode evolves towards 
a stationary mode of the lowest branch [see Fig. [l7|(b)]. 
In this case a substantial amount of power is radiated 
away due to the gap between the two branches. 



B. Mixed impurity that supports linear modes 

When both the bulk and defect have a self-defocusing 
nonlinear response {p = —1, f3 < 0) low-intensity lo- 
calized modes can still exist when the linear attraction 
(a > 0) dominates the nonlinear de- localization effect. 
In this case the nonlinearity-induced frequency shift is 
negative (see Sec. IVC) and the localized solutions exist 
for frequencies below the linear cut-off, uj < loq. The cor- 
responding dependence P{uj) is single- valued, as demon- 
strated in Fig. Hq. In the limit w — > 0, the power remains 
finite if cr < 2, but is unbounded otherwise. It follows 
from Eq. (|8|) that d < Go < Gf, and therefore ac- 
cording to Table || the corresponding localized modes are 
always stable. 



o 




FIG. 18. Power vs. frequency diagram for cr = 1, 7 = 1, 
a = 2, and (3 = —1. Dotted curve; asymptotic given by the 
power Pi{(^) of linear impurity modes. 



C. Mixed impurity Avith attractive nonlinearity 

Finally, we consider the case when localized modes do 
not exist in the linear limit, i.e. when a < and /3 > 0. 
As was pointed out above, the properties of localized 
modes depend on the relative strength of the bulk and de- 
fect nonlincarities, characterized by the ratio F. We have 
analyzed Eq. ( |27| ) and found that in the case of a "strong" 
defect with F > 1/2 the modes exist for all frequencies 
w > 0, and correspond to a single branch in the P(w) 
diagram, which asymptotically approaches Pi{(^) at high 
frequencies, according to Eq. ( p^ ) [see Fig. |l9|(a)]. Thus, 
these modes are stable above a certain cut-off, oj > cug, 
only if the impurity supports stable modes, i.e. for 7 < 1. 



This parameter ran ge c orresponds to region I in the dia- 
gram shown in Fig. ^. Note that for uj ^ ujq the linear 
impurity response is negligible and the mode character- 
istics should not depend on the sign of a. Indeed, we 
see that region I is the same in the cases shown in both 
Fig. |l| and Fig. 1^. 
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FIG. 19. Power vs. frequency diagram for defocusing bulk 
media with a — —2 and P — 3. (a) ct = 1, 7 = 0.8; (b) cr = 3, 
7 = 1, /3o — 2.38. Dotted curve: asymptotic given by the 
power Pi{io) of linear impurity modes. 
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FIG. 20. Stability regions for localized modes with a < 
and /3 > in a defocusing medium: I: the modes are stable 
for frequencies above a certain threshold; II: localized modes 
can exist and be stable in a bounded region of frequencies. 
In the region without shading localized modes exist for all 
u! > 0, but are unstable. The points (a) and (b) correspond 
to Figs. [19(a) and [l9|(b), respectively. 

For F < 1/2 (region II in Fig. pOl) the modes only exist 
if the nonlinearity of the impurity exceeds the threshold 
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value defined by Eq. (p6|), i.e. for (3 > /Sq. Then P{uj) has 
two branches, which appear a,i uj — Q and merge again 
at the critical point Ws, as shown in Fig. f9^(b). From 



III 



the properties of the critical points discussed in Sec 
we conclude that the stable modes correspond to (i) the 
upper branch (0 < a; < ujs) and (ii) part of the lower 
branch close to tOg, for which the slope dP/duj is positive. 



VII. COLLAPSE DYNAMICS 



A. Virial relation 



The nonlinearity-induced energy localization is a fun- 
damental physical problem. The localization can occur in 
the form of stationary nonlinear impurity modes, and it 
is essential to understand the underlying physical mech- 
anisms leading to localized states. A very efficient way 
for energy localization is the so-called collapse (or ''blow- 
up") dynamics when, under certain conditions, nonlinear 
self-focusing dominates over diffraction, leading to an un- 
limited growth of the field intensity in a finite time. In 
real physical systems, the actual "blow-up" never occurs, 
however, the initial collapse dynamics can be correctly 
described in the framework of the continuous model equa- 
tions, as long as the corresponding assumptions arc not 
violated. 

In this section we briefly discuss the collapse dynamics 
in the frame of our model equations ffl),(0) and derive 
sufficient conditions for collapse in the presence of an 
impurity. For the sake of simplicity, we are considering 
the case of the power-law nonlinearities, as introduced in 
Sec. WU- 



Studying the collapse-induced effects due to a nonlin- 
ear impurity, we distinguish two cases: (i) collapse in a 
bulk medium, away from the impurity site, and (ii) col- 
lapse at the impurity site. In the first case, the impurity 
acts similar to a small perturbation, and thus collapse can 
occur if the power of the nonlinearity exceeds the criti- 
cal value for a homogeneous self-focusing (i.e. p — +1) 
bulk medium, CTcr = 2 ||2^. We note that, for the case of a 
nonlinearly repulsive impurity (e.g., /3 < 0) at high inten- 
sities, the maximum field amplitude is always achieved 
away from the impurity site so that collapse can only 
occur in the bulk, see, e.g.. Fig. ||(b). Collapse at the 
impurity site can only take place if the impurity pos- 
sesses attractive nonlinearity (i.e. /3 > 0). We hereafter 
consider this case. 

In order to analyze the collapse conditions, we assume 
that the initial profile of the impurity mode is symmet- 
ric, and therefore -Hx) — ip{—x) at alH > due to the 
symmetry of Eq. (|l|). Then, the effective mode width R 
can be defined as follows: 



R^{t) = 



+ OC 



x^\^{x,t)\^dx. 



(37) 



We determine the temporal evolution of R{t) by follow- 
ing the standard procedure |2^ and derive the so-called 
virial relation: 



P 



d^R^ 



Ap 



dt^ 
(2 -a) 



= 8H + AaIo+ 



(1 + ^)7- 



+ 00 



m 



2cr+2 



da;-)- 4/3 



(1-7) 
(1+7) 



(38) 



7-7+1 



where P and H are the power and Hamiltonian defined 
by Eqs. (|) and (|), respectively, while Iq = |V'(0,i)P 
here designates the intensity at the impurity site. Be- 
cause P and H are conserved quantities, Eq. ( |38| ) can be 
integrated when the powers of the nonlinearities attain 
the critical values: ctci — 2, 7cr = 1, and linear defect 
response vanishes (a = 0): 



R{t) = jRl + AtHH/P), 



(39) 



where Rq is the width of the input mode, and we assume 
that the initial mode profile does not have a phase mod- 
ulation. If iy < 0, it immediately follows that the mode 
width decreases and eventually vanishes at a finite time, 
indicating that the peak intensity goes to infinity since 
the power P is conserved, thus the energy collapses to 
a single point — the impurity site. Therefore, a nega- 
tive value of the Hamiltonian is a sufficient condition for 
collapse in this case. 

If the power of the impurity nonlinearity is enhanced, 
i.e. 7 > 7cr, the corresponding term on the right-hand 
side of Eq. (^8|) becomes negative, indicating an increase 
in the collapse growth rate. Increasing the power of bulk 
nonlinearity above the critical value, a > a^ for p = -Hi, 
or decreasing the self-defocusing response, a < a^ for 
p ~ —1, we observe a similar effect. In all these cases, 
the relation H < remains a sufficient condition for col- 
lapse. The same argument holds for a defect with linearly 
repulsive response, a < 0. On the other hand, as we have 
demonstrated in the previous sections, a linearly attrac- 
tive defect can support stable localized modes which do 
not exhibit collapse. In this case, and more generally 
when at least one of the last three terms on the right-hand 
side of Eq. ( p8| ) becomes positive, we can not directly use 
the virial relation to predict the collapse conditions. 



B. Collapse conditions 

In order to predict whether collapse is possible for ar- 
bitrary values of the nonlinearity parameters we may em- 
ploy the connection between the collapse dynamics and 
the properties of stationary localized modes. First, we 
note that for the stationary solutions the right-hand side 
of Eq. ( pSf ) is identical zero. For the critical case a pertur- 
bation of the solution that makes H negative will result 
in a collapse. Thus, the collapse may occur when the 
high-frequency symmetric localized modes are unstable. 
Indeed, in the limit of high intensities we can neglect the 
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linear impurity response by taking a = 0, and then an 
exact collapsing solution can be obtained (at a = 2 and 
7 = 1) in the form: 



0{x,t) = ujCX{t) - x^\{t)/{'iC), 



(40) 



where X{t) = C{t — io)""'^, G is an arbitrary positive 
constant, Iq is the collapse time, and u{x) is the pro- 
file of a stationary localized mode with the frequency 
Lo. Although the solution (40) is unstable and there- 
fore does not describe the actual collapse dynamics, it 
demonstrates a link between the stationary modes and 
the collapse phenomenon. Moreover, for a homogeneous 
self-focusing medium it was proved that an unstable soli- 
ton will collapse if its power is slightly increased (see, 
e.g., Ref. 129|). Then, we arrive at the following state- 
ment: Collapse at the nonlinearly attractive impurity can 
be observed if and only if the stationary symmetric impu- 
rity modes exist in the limit of high frequencies and their 
power is bounded from above, i.e. P{uj — > -j-cx)) < Pmax, 
where Pmax is a constant. We note, however, that the 
stationary modes themselves are not necessarily unstable 
in the case with the critical power of the impurity non- 
linearity, as their power may asymptotically approach a 
constant from below (see the example in Fi g. p|). Finally, 



using the mode properties derived in Sec. IV C, we de- 
termine the parameter regions where collapse can occur, 
and summarize these results in Table II] . 



TABLE III. Collapse conditions 


Impurity response 


Bulk nonlinearity 


Collapse 


arbitrary 


p — +1 and a >2 


bulk 


/3 > and 7 > 1 
P > and 7 < 1 


p — —1 and o" < 27 
p = -1-1 and a > 2 


impurity 


/3 < and/or 7 < 1 
/3>0 


p = — 1 and/or a < 2 
p = — 1 and (T > 27 


does not occur 



waveguide embedded in a bulk Kerr medium. When the 
field localization becomes very high, the delta-function 
can no longer be used to approximate the layer response, 
and then the nonlinear response should be modified as 
follows 



T{I;x) 



/, |a;|^> d/2, 

a + 131, \x\ < d/2, 



(41) 



where d is the layer width, and the layer response pa- 
rameters are related to those used in the delta-function 
approximation, ad = a and fid = f3. As can be seen from 
Fig. El^(b), a decrease of the mode width and a growth of 
the peak amplitude stop when the energy gets localized 
inside the nonlinear layer. 





C. Collapse suppression 

The infinite growth of the field amplitude that re- 
sults from the collapse dynamics would not occur in a 
real physical system. However, in the framework of the 
model (|l|),(|2|) we observe that collapse can be induced by 
an impurity even in a standard nonlinear Kerr medium, 
i.e., when cr = 7 = 1, as demonstrated in Fig. pl|(a). 
The final stage of self-focusing, preceding the "blow-up" , 
should be described by modified equations that take into 
account some additional physical effects which can no 
longer be neglected when the collapse is approached. Be- 
low, we consider two possible mechanisms for such a col- 
lapse suppression. 

The first example is self-focusing of spatial nonlin- 
ear guided optical modes supported by a thin nonlinear 




FIG. 21. Comparison of different types of the col- 

lapse-induced dynamics, (a) Impurity-induced collapse — 
the continuum model (hi) for a self-focusing Kerr medium 
{a = J — 1, a — 2, f3 — 0.55); (b) Collapse suppression for 
a layer of a finite width (d — 0.02); (c) Collapse suppression 
due to the model discreteness {h ~ 0.02). 



16 



As the second example, we study the energy locahza- 
tion in an intrinsicaUy discrete system such as a waveg- 
uide array. It has been demonstrated that in homoge- 
neous lattices the collapse and the infinite growth of the 
peak amplitude is always suppressed |Q. This happens 
due to the presence of a minimal transverse scale, the 
characteristic distance h between the lattice sites or, for 
the present example, the width of the individual waveg- 
uides in the array. Therefore, when the width of the lo- 
calized mode becomes comparable with h, the evolution 
can no longer be described in the framework of the con- 
tinuum approximation. Therefore, the original model (H) 
should be modified to take the form of the discrete NLS- 
type equation: 



.di>n , (V'n+1 + 4'n-l " 2V^„) 



at 



2/l2 



+ .F(/„;n)V'„=0, (42) 



where n is the site number. An example of the energy lo- 
calization in a discrete system is presented in Fig. |2l|(c). 
When the mode becomes narrow, its confinement is de- 
fined by the discreteness, and the final state corresponds 
to breather-like oscillations near a stable stationary so- 
lution of a discrete model. 

Finally, we note that the initial stage of the collapse- 
induced dynamics is the same for all three different mod- 
els, as is clearly seen in Fig. |2^. It is also important to 
mention that the power of the high-frequency station- 
ary localized solutions of the modified models (El]),(E^) 
is no longer limited, and thus the observation of collapse 
suppression compl etely a grees with our general criterion 
introduced in Sec. VII B, 
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FIG. 22. Evolution of the field intensity at the nonlinear 
impurity corresponding to the plots shown in Fig. |2l|(a) — 
solid. Fig. pi|(b) — dotted, and Fig. El|(c) — dashed, respec- 
tively. 



VIII. CONCLUSIONS 

We have analyzed spatially localized nonlinear modes 
supported by a point-like impurity, in the framework of 



the generalized nonlinear Schrodinger equation. We have 
considered three possible types of such nonlinear impu- 
rity modes, i.e. symmetric one- and two-hump modes and 
asymmetric one-hump mode, and described their regions 
of existence and stability, for both focusing and defocus- 
ing nonlinearity of a bulk medium and two different (at- 
tractive or repulsive) types of the impurity. In particular, 
we have obtained an analytical stability criterion for non- 
linear localized modes based on the results of the linear 
stability analysis of the generalized NLS equation. For 
more specific physical applications, we have presented a 
detailed analysis of the nonlinear impurity modes and 
their stability in the case of the power-law nonlinearities 
in both the medium and defect, and discussed several 
scenarios of the instability-induced dynamics of the non- 
linear impurity modes. In particular, we have described a 
novel physical mechanism of the energy localization due 
to the impurity-induced collapse of a nonlinear mode at 
the defect site, which can occur when the power of non- 
linearity in the defect exceeds a critical value (i.e., 7 > 1); 
this effect can be observed for the Kerr medium as well. 
The problem we have analyzed above has a number of 
important physical applications ranging from the nonlin- 
ear dynamics of solids to the theory of nonlinear pho- 
tonic crystals and waveguide arrays in nonlinear optics. 
In particular, our results can be linked to different special 
cases of the theory of nonlinear guided waves in layered 
dielectric media, and they provide also a generalization 
of the theory of nonlinear impurity modes in solids, to- 
gether with a systematic classification of nonlinear im- 
purity modes and the analysis of the mode stability and 
its instability-induced dynamics. Additionally, this prob- 
lem can be considered as one of the first steps towards 
a deeper understanding of the competition between two 
different physical mechanisms of energy localization. In 
particular, we have presented our results emphasizing the 
cases where one can observe a clear evidence of compe- 
tition between the disorder- and nonlinearity-induced lo- 
calization. 
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